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ABSTRACT 



Context. Polarization of maser emission contains unique information on the magnetic field in the densest regions of massive star 
formation. 

Aims. Recently, the magnetic field induced Zeeman splitting was measured for the strongest known 6.7 GHz methanol maser, which 
arises in the massive star forming region G09.62+0.20. This maser is one of a handful of periodically flaring methanol masers. 
Magnetic field measurements can possibly provide insights into the elusive mechanism responsible for this periodicity. 
Methods. The 100-m Effelsberg telescope was used to monitor the 6.7 GHz methanol masers of G09. 62+0.20, in weekly intervals, 
for just over a two month period during which one of the maser flares occurred. 

Results. With the exception of a two week period during the peak of the maser flare, we measure a constant magnetic field of 
S|| » 1 1 + 2 mG in the two strongest maser components of G09. 62+0.20 that are separated by over 200 AU. In the two week period 
that coincides exactly with the peak of the maser flare of the strongest maser feature, we measure a sharp decrease and possible 
reversal of the Zeeman splitting. 

Conclusions. While the two phenomena are clearly related, the Zeeman splitting decrease only occurs near the flare maximum. 
Intrinsic magnetic field variability is thus unlikely to be the reason for the maser variability. The exact cause of both variabilities is 
still unclear, but it could be related to either background amplification of polarized emission or the presence of a massive protostar 
with a close-by companion. However, the variability of the splitting between the right- and left-circular polarizations could also be 
caused by non-Zeeman effects related to the radiative transfer of polarized maser emission. In that case we can put limits on the 
magnetic field orientation and the maser saturation level. 
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1. Introduction 

Polarization observations of astrophysical masers provide im- 
portant insights into the magnetic field properties of, among 
others, the dense r egions surrounding massive protostars (e.g. 
IVlemmingsi I2007L and references therein). In these regions, 
magnetic fields may play an crucial role in, e.g., suppressing 
fragmentation, altering feedback processes and stabilizing ac- 
cretion disks. Linear polarization observations of maser emis- 
sion can reveal the magnetic field morphology, while observa- 
tions of the circular polarization, generated due to Zeeman split- 
ting, can be used to measure the line-of-sight magnetic field 
strength. Especially for the non-paramagnetic molecules such 
as SiO, H2O and methanol, maser polarization also depends on 
intrinsic maser properties that determine the maser saturation 
level, such as brightness temperature, beaming angle and the rate 
of maser stimulated emission. Thus, polarization observations 
can, in addition to the magnetic field strength and structure, also 
provide constraints on mas er propert i es tha t are otherwise hard 
to determine (e.g. Vlemm ings et al.L l2006h . Circular polariza- 
tion, or Zeeman splitting observ ations have been mostly focused 



on OH and H2O masers (e.g. Hutawarakorn & Co henl 1999 : 



Sarma et al.1 120011; iBartkiewkz et all 120051; IVlemmin gs et al 



12006 !). However, recent observations have revealed significant 
Zeeman splitting of the 6.7 OHz 5i - 6oA^ methanol transition 
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(IVlemmingsl l2008', hereafter V08) in a sample of 17 out of 24 of 
the brightest northern methanol maser sources, indicating an av- 
erage magnetic field strength in the maser region of |B| ~ 20 mO. 

The massive star forming region O09.62+0.20 harbors the 
strongest known maser at 6.7 OHz. This maser, and its 12.2 OHz 
coun terpart, undergoes per i odic fl ares with a period of 244 days 
(e. g. lOoedhart et all I2003L l2004 . At the height of its flare, the 
6.7 OHz maser has been seen to reach a peak flux density of over 
7000 Jy beam"^ There is a time-delay of ~25 days between the 
flare in the strongest 6.7 OHz feature (hereafter the main fea- 
ture) at ~ 1.2 km s^'and the secondary feature at ~ -0.1 km s"', 
which also displays a different flare profile. The origin of the 
periodic behavior, however, is still unclear. The O09.62+0.20 
star forming region consists of a complex of HII regions in 
various evolutionary stages. Its 6. 7 OHz methanol m asers have 
been mapped with the ATCA by iPhillips et all (11998.) . and the 
strongest features are shown to b e associated with t he hyper- 
compact HII region labeled E by lOaray etaP (11993"). This re- 
gion is speculated to be excited by a BO star (Hofner et aI.L 
[l996h . VLBA obs ervations of the 12.2 OHz masers during the 
course of a flare ( Ooedhart et al.L l2005 l) indicate that the maser 
regions simply brighten in intensity, with no change in morphol- 
ogy, implying that the cause of the flare arises beyond the maser 
region. It has been previously speculated that the periodicity 
could arise from either the background HII region, or an infrared 
pump source. The monitoring observations presented here were 
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Fig. 1. The 6.7 GHz methanol maser total intensity (I) spectram 
of G09.62+0.20 for two observing epochs. The peak flux at July 
4 is scaled to the flux at Aug 14 2008 (top). The middle panel 
highlights the many weaker maser features of this region. The 
bottom panel shows the circular polarization spectrum (V) for 
the observations at Aug 14 2008. The thick solid line is the best 
fit fractional total intensity derivative. 

prompted by the detection of a possible magnetic field reversal 
between the main (By -3 mG) and secondary (By ^ 9 mG) 
maser features in V08. 

2. Observations and error analysis 

The 6668.519 MHz (5i -6oA^) methanol maser line of the mas- 
sive star forming region G09.62+0.20 was monitored weekly 
over a period of slightly more than 2 months (11 sessions from 
2008 June 14 to August 20) using the 5 cm primary focus re- 
ceiver of the 100-m Effelsber^ telescope. The observing dates 
were chosen to encompass the expected flare events of the main 
and secondary variable maser features. Unfortunately a setup 
problem made two sessions (June 14 and July 31) unusable. In 
addition to G09. 62+0.20, we observed G23.01-0.41 as a con- 
sistency check. As in the previous Effelsberg methanol maser 
Zeeman splitting observations (V08), the data were taken in 
position switch mode with a 2 minute cycle time. Data were 
collected using the fast Fourier transform spectrometer using 
two spectral windows, corresponding to the right- and left- 
circular polarizations (RCP and LCP). The spectral windows of 
20 MHz were divided in 16384 spectral channels, resulting in a 
~ 0.055 km s 'channel spacing and were centered on the local 
standard of rest (LSR) source velocities. The data were reduced 
as described in V08, with amplitude calibration performed on 
3C286. The spectrum of G09.62+0.20 for two Eff^elsberg obser- 
vational epochs is shown in Fig.[T] along with the circular polar- 
ization spectrum for one of the epohcs. As seen in the figure, the 
shape of each maser feature is identical at each epoch. However, 
the velocity of the maser peak of the main feature shifts similarly 

' The 100-m telescope at Effelsberg is operated by the Max-Planck- 
Institut fiir Radioastronomie (MPIfR) on behalf of the Max-Planck- 
Gesellschaft (MPG) 
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Fig. 2. The ratio A of the true and predicted channel rms noise 
as a function of channel flux for two epochs. The solid line indi- 
cates a two-component description of the rms noise increase due 
to dynamic range limits of the spectrometer. The dashed line in- 
dicates the anomalous increased rms noise for the maser features 
at Vlsr = 5.4 km s^'seen in the first epoch. 



in RCP and LCP by up to ~ 10 m s This is likely caused by 
a delay of the flare of weaker maser features in the wings of the 
main feature (Goedhart et al., in prep.). 

The target sources G09.62h-0.20 and G23.01-0.41 were ob- 
served for 10 and 6 minutes respectively. Due to changing con- 
ditions, the rms noise levels varied over the different monitor- 
ing sessions. Additionally, the LCP rms noise level was be- 
tween 50% and 300% larger than the RCP rms noise level. The 
LCP rms noise, which varied for the 10 min observations of 
G09.62-(-0.20 between 100 and 350 mJy, was thus the Hmiting 
factor to our Zeeman splitting determinations. Unfortunately, as 
the rms noise level was a factor of ~ 2 increased over that in 
the V08 observations, we were only able to detect significant 
Zeeman splitting in our consistency check source G23.01-0.41 
at 4 of the epochs. Within the errors, these epochs were consis- 
tent with the V08 observations. 

Synchronous observations of G09.62+0.20, as part of a 
program to monitor the variable methanol maser sources 
(Goedhart et al., 2005), were carried out using the 26m telescope 
at Hartebeesthoek Radio Astronomy Observatory (HartRAO). 
Full sampling of the flare was not possible due to schedul- 
ing constraints. Flux calibrations were done using continuum 
drift scans across Hydra A and 3C123 and a further check was 
made using the methanol maser G351.42-h0.64 as a comparison 
source. Comparison with the HartRAO measurements enables 
us to estimate the absolute flux errors to be less than 10%. 

As our observations were aimed at detecting possible small 
variabilities of the Zeeman splitting, we performed an addi- 
tional analysis of the channel rms noise in the spectrum of 
G09. 62-1-0. 20. To determine the increase of channel rms noise 
as a function of maser flux, we used the 5 individual 2 minute 
scans of two of the epochs to determine the channel rms in those 
channels that contained significant maser emission. The individ- 
ual maser spectra were normalized to the peak flux of the com- 
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bined scans to minimize the effect of intrinsic maser variability. 
We then calculate Ai = rmsi/cr for each channel /, where cr is the 
rms noise value for the emission free channels. Fig.|2]shows A 
as a function of the maser flux in each of the channels with > 5cr 
maser emission. We find that for all of the maser features, with 
one exception, A stays approximately constant within a factor of 
~ 3 up to a maser flux of ~ 50 Jy beam"' after which it increases 
with A oc (Flux)" ''. During the first epoch, only the maser fea- 
ture at Vlsr = 5.4 km s^'did not follow this relation and already 
deviated from the expected rms noise level when its flux became 
> 20 Jy beam As this did not occur during the last epoch, this 
is possibly due to weak narrowband interference. Since the noise 
characteristic is similar for both polarizations and all epochs, the 
rms error increase is unlikely to be due to receiver saturation 
and would be unable to cause a systematic shift between RCP 
and LCP. However, the analysis does imply that the errors on 
Zeeman splitting determination in V08 should be increased. For 
the majority of the masers in V08, this increase is less than a fac- 
tor of ~ 5. However, for the masers of G09.62h-0.20 the errors 
need to be increased with a factor of 40 and 8 for the main and 
secondary maser features respectively. 



3. Variability of RCP-LCP frequency splitting 

3.1. The case of G09.62+0.20 

We determined the Zeeman splitting of the two brightest 
6.7 GHz methanol maser features of G09.62h-0.20 using the 
RCP-LCP cross-correlation method as described in V08. Fig. [3] 
shows the line-of-sight magnetic field strength By for the two 
maser features at each epoch as derived from the Zeeman split- 
ting AVz, using 0.049 km s ' G ' as the best estimate for 
the Zeeman splitting coefficient. In addition to the monitor- 
ing epochs, the figure also includes the folded in observations 
presented in V08. We find that the line-of-sight magnetic field 
strength in the secondary maser feature is stable during the ob- 
servations, with an error weighted average magnetic field of 
B|l - 10.9 + 2.3 mG. In contrast, the Zeeman splitting of the 
main maser feature sharply decreases at the time the maser flare 
reaches peak flux. Monitoring at weekly intervals has proven to 
be too coarse to put strong constraints on the length of the pe- 
riod with decreased AVz- However, while the typical duration 
of a flare is up to two months, the decrease in Zeeman splitting 
lasts only for an approximate two week period around the peak 
of the flare. Determining a error weighted average magnetic field 
for the main maser feature using the 7 epochs on either side of 
the two week period with decreased Zeeman splitting, we find 
B|l = 11.0 + 2.2 mG. Thus, the magnetic field strength is re- 
markably similar on both maser s, even though both fea tures are 
separated by more than 200 AU jGoedhart et al.Ll2005l) . 

Fig. [3] thus indicates that, when the flare of the main maser 
feature reaches its peak flux, the Zeeman splitting decreases sig- 
nificantly and potentially even changes sign. As the observations 
of V08 that also revealed a much lower (and possibly reversed) 
AVz were taken close to the peak of the previously flaring pe- 
riod, this behavior appears to repeat itself regularly. No signifi- 
cant Zeeman splitting decrease is seen for the secondary maser 
feature. The observations with the HartRAO telescope showed 
that the flare of this feature, while typically approximately 25 
days after the flare of the main feature, was much more irregu- 
lar and reached its peak at approximately the same time as the 
main feature. However, throughout the Effelsberg observations 
its flux variations were only ~ 10% while that of the primary 
feature was over 20%. 
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Fig. 3. (bottom panels) The Zeeman splitting AVz (in m s"') and 
derived line-of-sight magnetic field strength By (in mG) for the 
two strongest 6.7 GHz methanol maser features of G09.62+0.20 
at the 9 successful monitoring epochs (filled dots). The previ- 
ous observations of Nov 1 1th 2007 (open square; V08) has been 
folded into the new observations. The vertical short dashed lines 
indicate the predicted date of the emission peak. The horizon- 
tal dashed lines indicate the weighted average magnetic field 
strength (see text), (top panels) HartRAO telescope observations 
of the two maser features for two flaring periods. The period 
of the Eff'elsberg observations is indicated by the vertical short 
dashed lines. 



3.2. Possible origin of tine observed variabiiity 

The cause for the sudden decrease in Zeeman splitting during the 
maser flare is unknown. We can confidently rule out instrumental 
effects as the reason for the decrease of AVz of the main maser 
feature for a number of reasons. First of all, no significant corre- 
sponding decrease is found for the secondary maser feature at the 
same epochs. Furthermore, the quite constant Zeeman splitting 
measured before and after the peak, and the fact that the negative 
Zeeman splitting was also found during the maser flare 8 months 
earlier, point to the stability of the instrumental setup as well as 
the robustness of the data reduction and analysis method. 

The observed eff'ect is thus likely intrinsic to the source. The 
measured decrease in AVz however, starts after the maser flux 
has already entered the flaring stage. Thus, while the measured 
Zeeman splitting variation is related to the maser flare, it is un- 
likely that the flare itself is caused by changes in the magnetic 
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field. Here we describe three possible scenarios that could give 
rise to the observed effect. 

Background amplification: As the observed Zeeman split- 
ting is generated by the average magnetic field throughout 
the entire maser region, the most straightforward cause of a 
drop in observed Zeeman splitting is the superposition of two 
maser regions wit h opp osite magnetic fields. As described in 
iBoboltz et al.1 (e.g. lT998 ). maser flares can often be attributed to 
the chance alignment of maser regions, when a foreground maser 
amplifies the emission from a background maser Alternatively, 
instead of a maser, the background source could be a strongly 
polarized continuum source. In this case, polarization of the 
seed emission can cancel out any circular polarization gener- 
ated within the maser. However, while this simple model also 
naturally explains the maser flare, the origin of the observed pe- 
riodicity remains unclear 

Intrinsic magnetic field variability: A second option is that 
the observations reveal an actual change in the magnetic field 
within the maser region. Various mechanisms can produce a 
change in magnetic field. It has been suggested that the peri- 
odic maser flares are due to changes in the maser pumping re- 
sulting from a binary interaction betw een a massive protostar 
and companion dGoedhart et al.Ll2b03h . Interactions between the 
magnetospheres of the two companions could cause periodic be- 
havior of the magnetic field, as has been observed in young low- 
mass binary systems (e.g. Massi et al., 2008). However, typical 
magnetic reconnection events do not last long enough to explain 
the approximate two weeks duration of the magnetic field vari- 
ability. It is also unclear if the magnetic interaction would be 
noticeable in the methanol maser region at several hundred AU 
distance from the protostars. Still, an embedded binary will also 
give rise to other complex interactions besides that of the mag- 
netic field, such as those between possible accretion disks and 
outflows. These interaction will be imprinted onto the observed 
magnetic field. Thus, if the observed magnetic field variability is 
truly due to intrinsic changes of the magnetic field, any possible 
explanation of the periodic maser flares will also need to take 
into account the behavior of the magnetic field. 

Maser radiative transfer: Alternatively, we are observing 
a combination of Zeeman and non-Zeeman effects, competing 
when the maser is at its brightest. Non-Zeeman effects were 
briefly discussed in V08. Specifically the effect where the axis 
of symmetry for the molecular quantum states rotates when the 
maser saturates, and the rate of stimulated emission R becomes 
larger than the Zeeman frequency shift gO., bears further inves- 
tigation. An unfortunate error was introduced in the calculation 
of gQ presented in V08, with the true gO. being smaller than pre- 
sented there. For a typical magnetic field strength B, in the dense 
maser region, of order 1 mG, gD. ~ 13s', approximately three 
times larger than the rate of stimulated emission R ~ 4 s ' for the 
most saturated 6.7 GHz methanol masers with a maser bright- 
ness temperature Tb ~ 10'^ K sr. This assumes a typical maser 
beaming angle AO = 10"^ sr, which decreases rapidly when the 
maser saturates. Thus, for most typical methanol masers with Tb 
of order a few times 10'" K, R < 0.1 s ' << gQ. and little or no 
intensity dependent polarization is generated. This is supported 
by the fact that no relation between maser intensity and Zeeman 
splitting was found in V08. However, the brightest masers of 
G09. 62-1-0. 20 could have an R that approaches or even becomes 
larger than gQ, raising the possibility of intensity dependent cir- 
cular polarization mimicking the Zeeman splitting between the 
RCP and LCP spectra. The generation of circular polarization 
due to t his effect has been inves tigated for a. J - 2 - \ transi- 
tion bv iNedoluha & Watsonl([T990i) . Although the effect is likely 



smaller for the transition of the 6.7 GHz methanol masers in- 
volving higher angular momentum states, it is found that the 
sign of the circular polarization, and consequently the sign of 
the splitting between RCP and LCP, is opposite from that gen- 
erated by the regular Zeeman effect when the angle between the 
magnetic field and the maser line of sight obeys sin^ 6 < 2/3. 
Observationally, this implies that with increasing maser bright- 
ness temperature and consequently R, the observed splitting be- 
tween RCP and LCP spectra decreases. As the average Zeeman 
splitting of the two maser features, with fluxes different by an 
order of magnitude, are identical, the non-Zeeman effect appar- 
ently only becomes important when the flare nears its maximum 
flux. This indicates the masers of G09. 62-1-0. 20 approach com- 
plete saturation at its peak. It will be possible to test the hypoth- 
esis of the non-Zeeman interpretation of the observed splitting 
variability by simultaneously observing the maser linear polar- 
ization, as similar considerations predict an intensity dependence 
of fractional linear polarization and polarization angle. 

4. Conclusions 

We have presented Effelsberg 100-m telescope monitoring ob- 
servations of the RCP-LCP frequency splitting of the periodi- 
cally flaring 6.7 GHz methanol masers of the massive star form- 
ing region G09.62+0.20. A significant decrease in Zeeman split- 
ting and thus possibly magnetic field strength is detected on the 
strongest maser feature during a two week period surrounding 
the flare maximum. Besides this decrease, a remarkable constant 
B\\ of ~ 1 1 mG is detected for the two brightest maser features 
separated by over 200 AU. The cause for the decrease in mea- 
sured Zeeman splitting is still unclear but it is either related to 
the mechanism that causes the periodic maser flaring or a result 
of a non-Zeeman effect when the maser saturation level becomes 
significant. 

Thus, G09. 62-1-9. 20, with the relative predictability of the 
flares, varying levels of intensities amongst maser features and 
correlated variability in RCP-LCP frequency splitting, is an ideal 
natural laboratory to test theories relating to both massive star 
formation and maser physics. 
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